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We propose a class of networks which can be regarded as an extension of the graphitic network. 
These networks are constructed so that surface states with non-bonding character (edge states) are 
formed in a tight-binding model with one orbital for each atomic site. Besides, for several networks, 
the tight-binding electronic structures become a zero-gap semiconductor. These properties have 
been found in the vr-electron system of the graphene. Thus, we call these networks hypergraphite. 



I. INTRODUCTION 

Graphite has many intriguing aspects from scientific 
and industrial view points. It has high electronic con- 
ductivity and large anisotropic diamagnetic susjceptibil- 
ity. The electronic structure is semi-metallicE'EI An im- 
portant point is that these characteristic features of an 
electronic structure originate from a particular lattice 
structure of graphite. The lattice structure is a stack 
of honeycomb networks of carbon atoms. A single layer 
of graphite is called graphene. 

Our interest in graphite has deepened much by findings 
of rich physics not only in pure graphite itself but arlso in 
various graphitic materials. For example, fuUerenesu and 
carbon nanotubesu are regarded as a deformed graphitic 
structure. The discoveries have revealed that the topol- 
ogy of sp^ carbon networks crucially influences their tt- 
electronic structures. An example showing effects of 
topology is that electronic structures of single -wall car- 
bon nanotubes depend on the chiral vectorsJj'B The tt- 
electronic structures show metallic or semi-conducting 
behavior. Interestingly, recent experiments confirmed 
this interplay betweenJie topology of sp^ carbons and 
the TT-electronic states Ela 

Let us start from the graphene to understand the elec- 
tronic structures of graphite and related materials. The 
electronic structure of the graphene around the Fermi en- 
ergy is composed of a tt (or 7r*)-band. The 7r-band and 
the 7r*-band are the highest valence band and the lowest 
conduction band, respectively. The two bands degener- 
ate at the edge of the first Brillouin zone (1st BZ), called 
K (or K')-point, and show /c-linear dispersion at these 
points. Thus, the graphene is a zero-gap semiconductor. 
Namely, the density of states (DOS) becomes zero at the 
Fermi energy and linearly increases with leaving from the 
Fermi energy. 

Another interesting 7r-electron system is nanometer- 
sized graphite called "nanographite" .Q Nanographites are 
a class of mesoscopic systems which are situated between 



aromatic molecules and graphite. In nanographite, the 
presence of edges, their shapes and size crucially affect 
their 7r-electronic structures. 

We introduced ribbon models of a graphene to study 
effects of the edge on their 7r-electronic structures theo- 
retically. The typical edge shapes are a zigzag edge or an 
armchair edge. Hence, we investigated zigzag ribbons, 
which are graphene ribbons with only zigzag edge, and 
armchair ribbons, which are graphene ribbons with only 
armchair edge. 

There is a clear difference between 7r-electronic struc- 
tures of these two ribbon models. In case of zigzag rib- 
bons, a localised, states appear around the edge at the 
Fermi energy.Elil^ While such localized states do not ex- 
ist on the edge of armchair ribbons. Thus appearance of 
the localized states is a topological effect of zigzag rib- 
bons. 

The localized states on zigzag ribbons were named 
"edges states" , because they exist around edges of zigzag 
ribbons. A solution of the edge states was constructed 
on a semi-infinite graphene with a zigzag edgca This 
solution shows that the edge states is a non-bonding 
orbital (NBO). Band structures of zigzag ribbons pos- 
sess a pair of almost fiat bands near the Fermi en- 
ergy. These flat bands, i.e. edge states, induce a sharp 
peak in DOS at the Fermi energyllj. Therefore, zigzag 
ribbons are expegted to show unusuaLjaroperties, e.g. 
spin poiarizationjj lattice defcManationJiJ magnetic field 
effects,t3 transport propertiesJl3 etc. The appearance of 
edge states on zigzag ribbons is also recognized in terms 
of the first-principles calculations within the framework 
of the local density approximation.EjH 

Now it is natural to ask a question whether other struc- 
tures which have edge states due to a topological reason 
exist or not. This is because interesting phenomena have 
been found only for zigzag ribbons of a graphene. 

In this paper, we report on a method to construct an 
extension of the graphitic network. The extended net- 
works possess the same characteristics that the graphene 
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shows both in the topological network and in the elec- 
tronic structures. The networks are an AB bipartite net- 
work. When a proper boundary condition is subjected to 
them, edge states appear in these networks. Besides, the 
bulk states of these networks become a zero-gap semi- 
conductor. And the highest valence band and the lowest 
conduction band degenerate at the Fermi energy in 1st 
BZ and the two bands show /c-linear dispersion at the 
Fermi energy. We name the networks "hypergraphite" , 
because they are regarded as an extension of the graphitic 
network. 

The organization of this paper is the following. In sec- 
tion II, we argue about edge states. In particular, we 
introduce networks in which edge states appear under 
a proper boundary condition. In section III, we argue 
about bulk states of these networks by using a numeri- 
cal calculation and an effective mass approximation. In 
section IV, we give a definition of the hypergraphite. 



II. EDGE STATES 

The purpose of this section is to show several networks 
in which edge states appear. First, we recall a solu- 
tion of edge states on a semi-infinite graphene with a 
zigzag edge. Then we consider the relationship between 
the topology of zigzag ribbons and the solution of edge 
states. This consideration leads us to find a method to 
construct these networks. Next, we give solutions of edge 
states which appear in these networks. We show band 
structures of slab model of them. 

In this section and the next section, we use a single- 
band tight-binding model (s-TBM). This is because we 
study the relation between topology of networks and elec- 
tronic structures. The Hamiltonian is written as 

where c] ^{ci^a) creates (annihilates) an electron with 
spin, fj, on the i-th site (atom). We let the on-site energy 
zero, which is always possible for s-TBM. A bond con- 
nection is represented by and taken only between 
connected sites with hopping integral, t. Hereafter we as- 
sume that the hopping integral, t, is unity for simplicity. 
And we consider only a paramagnetic state. In addition, 
each site is occupied by one electron on the average. 



A. EDGE STATES ON A ZIGZAG RIBBON 

We recall how to construct the solution of the edge 
states on a semi-infinite graphene with a zigzag edge.cl 
The solution, (pk,D{x,z), is for the 7r-electron systems 
(Fig. 1(a)). Here, k and D denote a wave vector along the 
edge (cc-direction) and a dumping factor inward the rib- 
bon (z-direction), respectively. Each atom at the zigzag 



edge is assumed to be terminated by a hydrogen atom. 
The honeycomb network is AB bipartite and each atomic 
site can be classified into A- or B-site. In terms of this 
classification, all of the zigzag-edge sites of a semi-infinite 
graphene is in one of the two sublattices. In this paper, 
we define sites in the sublattice containing the edge sites 
as A-sites and sites in the other sublattice as B-sites. 

We list some characteristics of (f>k,D{x,z). 1) The am- 
plitude is non-zero only on A-sites. This means that edge 
states are NBOs. 2) Edge states are a dumping wave with 
a dumping factor, D — 2 cos I, when ^ < fc < tt. 3) 
The edge state completely localizes at the edge, {D = 0), 
when k = TT. The amplitude are -1-1 and -1 alternately 
along the edge except for normalization. 4) The edge 
state approaches a propagating wave {D 1), when 
fc — !■ This limit, 0^2— >7r/3,D— >ii coincides with bulk 
states at the K point in BZ of the graphene. 

A construction method of the solution of the edge 
states in a semi-infinite graphene with a zigzag edgea is 
briefly summarized as follows. We use a notation (x, nj), 
{j = A OT B), to specify a site. Here, x denotes a posi- 
tion of a cell represented by dashed lines in x-direction 
(See Fig. 1(b)). The cell contains two types of sites of 
both A- and B-sites. nA (or nB) represents A sites (or B 
sites) in the n-th zigzag-chain in the cell, respectively. We 
start from setting a trial function, ipkix, lA) = exp{ikx), 
to construct the solution of edge states. Next, we de- 
termine amplitudes of ipki^, iA)^ (i = 2, • • • , oo), so that 
'ipkix,jB),{j = l,---,co), become zero. The final ex- 
pression becomes an eigen function with zero amplitude 
on B-sites, i.e. edge states. 



B. NETWORKS WITH EDGE STATES 



We consider the relationship between topology of 
zigzag ribbons and the above method to construct the 
wave functions. Usually, a network of zigzag ribbons 
is regarded as assembly of hexagons. However, we give 
another topological view point to the network of zigzag 
ribbons. It is that the network of zigzag ribbons are re- 
garded as a network made by linking A-site of the i-th 
AB bipartite zigzag-chain and B-site of the {i + l)-st one 
by extra bond (i = 1, • • • , oo) (Fig. 1(c)). A character- 
istic of the view point is to connect zigzag-chains which 
are an AB bipartite lattice and have a NBO as an eigen- 
function. 

We conjecture conditions for networks having edge 
states from the topological view point as follows. 1) A 
network is regarded as one made by successively con- 
necting units. The unit is an AB bipartite lattice which 
has a NBO as an eigen-function. 2) The connection is 
made by linking A-sites of the i-th bipartite lattice with 
B-sites of the (i + l)-st bipartite lattice by extra bonds 
(i = 1,2, • • -,00). 
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Fig.l. (a) A solution of an edge state for a semi- 
infinite graphene with a zigzag edge. Closed circles 
and opened circles represent A-sites and B-sites, re- 
spectively. The wave function has finite amplitude 
at A-sites. D is a dumping factor, D = 2cos^. 
(b) This figure shows the notation which we use in 
text. The arrow represents a unit vector. Dash line 
shows a unit cell of a semi-infinite graphene with 
a zigzag edge. N represents the numbers of zigzag 
chains, (c) A way to construct a zigzag ribbon from 
zigzag chains. Black bonds are extra bonds to con- 
nect zigzag chains 





Fig. 2. (a) The three-dimensional three-fold coordi- 
nated network (the 3-3 network). N represents the 
number of blocks shown in Fig. 2(c). Closed cir- 
cles and opened circles represent A-sites and B-sitcs, 
respectively, (b) ID zigzag chains aligned in the x- 
direction. (c) A structure made by connecting zigzag 
chains aligned in x-direction and y-direction. 

Here, we introduce three examples which satisfy the 
above conditions. They are the three-dimensional three- 
fold coordinated network (the 3-3 network), the diamond 
structure (DS) and the three-dimensional five-fold coor- 
dinated network (the 3-5 network). We show the 3-3 
network in Fig. 2(a), the DS in Fig. 3(a), and the 3-5 
network in Fig. 4. Of course, edge states appear in all 
of them under a proper boundary condition as shown in 
the next subsection. 

We show that they satisfy the above conditions. The 
3-3 network is regarded as a network connecting AB bi- 
partite lattices as shown in Fig. 2(b). The AB bipar- 
tite lattice is made by combining a set of zigzag chains 
aligned in a;-direction (Fig. 2(c)) and in ?/-direction, al- 
ternately. A NBO exists on the AB bipartite lattice as 
an eigen function. The DS and the 3-5 network are re- 
garded as networks of connected honeycomb lattices and 
of connected square lattices, respectively. Naturally, a 
honeycomb lattice and a square lattice are an AB bipar- 
tite lattice and have a NBO as an eigen-function. 
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Fig. 2. (d) A solution of the edge state for the semi-infinite 3-3 network with zigzag surface. The wave func- 
tion of the edge states has finite amphtude at the sites indicated by closed circles. Dashed lines represent 
the next layer to the layer shown by solid lines. And two arrows represent two-dimensional unit vectors. 
The wave function is given by \E'^„ = D'-~^ exp{i{kxn+kym)), = D''~^{— 2 cos kx)exp{i{kxn+kym)) 

and D = 4 cos kx cos ky. D is a dumping factor, (e) The 1st BZ of a slab model of the 3-3 network. Edge 
states emerge in the shadowed region, (f ) A band structure and DOS of the 3-3 network. Here, we adopt 
a slab model with N = 10. 



C. SOLUTIONS OF EDGE STATES 

We show solutions of edge states which appear in the 3- 
3 network and the DS under a boundary condition. The 
boundary condition is an open boundary condition with a 
surface. The surface is composed of A-sites being only in 
a connected AB bipartite lattice used as a unit to make 
the full network. We call the surface "zigzag surface". 
A dumping factor, D, and an area in which edge states 



appear in 1st BZ arc determined from the solution. 

We show solutions of edge states in Fig. 2(d) for the 
3-3 network and in Fig. 3(b) for the DS. In this paper, 
we do not show a solution of the 3-5 network, but a sim- 
ilar solution is constructed for it. Here, we note that 
s-TBM is used when solutions of edge states are con- 
structed. Namely, we assume that a hopping integral, t, 
is the same and unity between any connected sites in any 
direction. 



4 



(a) (b) 




r M K r DOS 



Fig. 3. (a) A diamond structure. N represents the number of honeycomb lattices, (b) A solution of the 
edge states for semi-infinite DS i,e, with (lll)-surface. In each panel, dashed lines represent a next layer 
to a layer represented by solid lines. The wave function of the edge states has the amplitude at the sites 
indicated by closed circles. And two arrows represent two-dimensional unit vectors. The wave function 
is given by ^'^j^^ = (-D)'"^ exp(z(/ca;n + kym)) and D = eK^^fe^+fea) _|_ gi(fe,-2fc„) _|_ gi{k^+ky) ^ D is a 
dumping factor, (c) The 1st BZ of a slab model of DS. Edge states emerge in the shadowed region, (d) 
A band structure and DOS of a slab model of DS. Here, we adopt a slab model with N — 20. 



D. BAND STRUCTURES OF SLAB MODELS We recognize appearance of edge states from these re- 

sults. 



The electronic band structures of a slab model are 
shown in Fig. 2(f) for the 3-3 network and in Fig. 
3(d) for the DS. They are obtained by solving the eigen- 
value equations of s-TBM for each structure numerically. 
Since, we consider the case that the electron number is 
the same as the number of atomic sites, the Fermi energy 
becomes just zero for these modelalj. 

A pair of almost flat bands exists at the Fermi energy 
in both band structures. There is a sharp peak at the 
Fermi energy in DOS for the 3-3 network and the DS. 



Edge states appear in the 3-3 network and the DS 
which satisfy the conditions described in subsection B. 
In this paper, we show only two examples, but we have 
found many networks which satisfy the conditions and 
confirmed that edge states appear in these networks. In 
addition, the dimension of the structure does not have 
to be three or two. It is possible to design a higher di- 
mensional structure with edge states. A dumping factor 
for these networks is obtained by using a transfer matrix 
method (See Appendix. A). 
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Fig. 4. A three-dimensional five-fold coordinated net- 
work. N represents the number of square lattices. 



III. BULK STATES 



In this section, we argue about bulk states of the net- 
works which satisfy the conditions described in the last 
section. Electronic structures of the 3-3 network and 
the DS with a zigzag surface are similar to that of the 
graphene. Here, we have a question of how the electronic 
structure of bulk states is. Bulk states of the 3-3 network 
and the DS are made use of in this section as examples. 
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Fig. 5. (a) A band structure and DOS of a bulk of 
the 3-3 network, (b) The 1st BZ of the 3-3 network. 
In this figure, a region which contain fc^ > and 
ky > is showed. The fc-line is represented by the 
thick line. 



A. BAND STRUCTURES 



We show electronic band structures and DOS in Fig. 
5(a) for the 3-3 network and in Fig. 6(a) for the DS, 
respectively. The Fermi energy is jiist zero in the model. 

At first, we explain the band structure of the 3-3 net- 
work. As shown in Fig. 5(a), the lowest conduction band 
and the highest valence band degenerate at a point on the 
T — X line. At this point, two bands have A;- linear disper- 
sion. This is the same characteristic that the 7r-states of 
the graphene have. But, in the case of the 3-3 network, 
the lowest conduction band and the highest valence band 
degenerate on a closed line in 1st BZ. This line is shown 
in Fig 5(b). Recall that the lowest conduction band and 
the highest valence band degenerate at only K (or K') 
point in the case of the graphene. This line for the 3-3 
network is understood to be an extension of the K point 
of the graphene. Hence, we call this line "fc-line" . This 
point is argued by using the k ■ p approximation in next 
subsection again. In addition, the system is a zero-gap 
semiconductor as seen in DOS. 

Next, we explain the band structure of the DS. There 
is a A;-line, where the lowest conduction band and the 
highest valence band degenerate and two bands possess 
fc-linear dispersion at W point. These fc-lines are shown 
in Fig. 6(b). The system is also a zero-gap semiconductor 
as seen in DOS. 



(a) 




L r X w DOS 




k - line 



Fig. 6. (a) A band structure and DOS of a bulk of 
the DS. (b) The 1st BZ of the DS. In this figure, a 
region which contain k^ > and ky > is showed. 
The fc-lines are represented by the thick lines. 
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The electronic structures of both the 3-3 network and 
the DS show the same characteristic. However, there is a 
difference in the shape of fc-hne of these two structures. 
The k-line is closed in 1st BZ in a case of the 3-3 network, 
while the fc-line is open in 1st BZ in a case of the DS. 

We have studied many networks which satisfy the con- 
ditions described in the last section except the 3-3 net- 



work and the DS. The most of the networks show the 
same characteristics that the 3-3 network or the DS 
shows. One point which we have to remark is that in 
some structures made by our construction method, the 
bulk band structure is not a zero-gap semiconductor but 
a metal. In these exceptional cases, edge states appear 
in these networks whose bulk states become metallic. 



B. K P APPROXIMATION 



In this subsection, we study bulk states of the 3-3 network and the DS using the k -p approximationc3. The results 
show again that there are bands with fc-linear dispersion in these networks. 

First, we study the 3-3 network by the k -p approximation. We start from obtaining wave functions on a fc-line by 
diagonalizing the single-band tight-binding Hamiltonian. Wave functions are represented in a Bloch form as 

^''('■^ = ^ E E *5^P(*k • R.)6,(k)</), (r - R,), (1) 



3 = 1 R; 

where (f>j(jc — R^) is a wave function of the j-th site in a unit cell. The vector, R.^, represents the position of the z-th 
cell. bjCk) is an amplitude which should be determined. iV is a number of unit cells. Here, we assume that 4>j{r — Ri) 
are orthonormal: 

(0j(r - Ri),0j/(r - Ri/)) = (5r,,r^,5jj'. 
The eigenvalue equation of the 3-3 network becomes following, 

i/(k)b(k) =eb(k). (2) 

The matrix, i?(k), is defined by, 

Hi^)=(l,o)^ (3) 



where 



js^^i -{I + cxp{i\^kxa) -exp{i{^kxa+ ^kya + Sk^a) i 

~^ -1 ^(1 + cxp(iV^fcya)) '■ ^' 



a is the bond length. A vector b(k) is given by 

b(k) 



/6i(k)\ 

&2(k) 
&3(k) 

V 64 (k) / 



We can deduce equations determining a fc-line on which energy is equal to zero from the above eigenvalue equation 
(Eq. (|)). The condition, detA = 0, gives fc-line equations, 

fc, = ' ' ' 

where K = {K^, Ky, 0) is a point on fc-line. 

We need wave functions on fc-line to apply the fc • p approximation. These are given in a Bloch form as, 

V'A(r) = -^^{exp(iK.RO(-expH/3))y(/)i(r-R,)+X02(r-R.))}, 
R. 

V'B(r) - ^= ^{exp(zK • K,)iYMr - " exp(m)X,/)4(r - R,))} , 
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where 



X — \J2 cos a, Y ~ ^2 cos/3, L = ^2 cos a + 2 cos/?, 

and 

a = —K^a, fJ = —Kya. 

We introduce a wave vector k measured from a point on fc-line and define two functions, 

$j-k(r) =exp(ikT)V'j(r), (6) 
with j — A, B. These functions are orthonormal: 

(*j,k,$j',k') = %"5kk'. 
The wave function near the point can be expanded as 



j=A.,B'' ^ ' J=A,B 

where Fj (r) is the envelop function defined by 

F^r)^ J ^exp(zk.r)Q(k). (8) 



Substituting Eq.(^ into the Schrodinger equation, we obtain the following k ■ p equation, 

f -^exp(i(a + /?))(- tan a 4 - ^atan/3 ky + i 3a k^)FA{v) = EFb{v) 
\ -^exp(-z(Q! + /3))(-^atanQ; 4 - ^atan/3 fey - * 3a 4)FB(r) = £;i^A(r) ' 



(9) 



where kx^y^z are defined as kx.y,z = '—i^x,y.z- Here, we rewrite the above expression with kj_ and fc||, which are 
components of the operator, k, perpendicular to fc-line and parallel to fc-line in k^ky-plane, respectively. They are 
given by 

k± \ _ f cos 6k sin 0k \ / \ / = 3^ tan a 

kii ) ^ \ -sinOK cosOk ) \ L J ' 1 sin^K = ^tan/? ' 



Ak = Y tan'' a + tan^ /3. 
By using the components, fcj^ and fcii , Eq. (W becomes 



exp(i(a + PMKk^ - iijkz) \ ( FaIv) \ ^^^'"^ ) (10) 



exp{~i{a + (3)){CKkj_ + irikz) ) \Fb{y) ) \Fb{v) 

where 

eK = -^aAK/i2 
I] = Za/F^ 

The final expression is the same as that of the k ■ p equation for the graphene around K (or K') point. We can obtain 
the energy eigenvalue around fc-line from which is 



E^^^^ki+rfkl (11) 

Here, k± and kz are a component of the wave vector, k, defined above. Thus, the energy dispersion around a fc-line 
is in proportion to fc, which is a distance from the fc-point to the fc-line in plane perpendicular to k-line. 
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Next, we study the DS by using the k ■ p approximation. We have several fc-hnes in 1st BZ in this case. We choose 
a fc-hne, 



Ky=0 

^ a 

to demonstrate the k ■ p approximation. From the same procedure, we can obtain the following k ■ p equation, 



(12) 



eM^^MKjy + iVK^k,) \ ( Fa{v) \_j,f FA{r) , , 



where 



= asm(^— ), r]K^ = acos(^— ). 

This equation is again the same form as that of the graphene. We can obtain the energy eigenvalue around the 
fc-line from Eq.(O), which becomes 



E=^eKkl+vlkl. (14) 

Since we consider k-line which is parallel to fc^-axis, this energy eigenvalue is proportional to a component k of the 
wave number, k, in a plane perpendicular to k-line, i.e the fcj,fc^-plane. The point (0,0, ^), where Eq.(^3|) possesses 
a singularity, is an exception. The energy eigenvalue is proportional to only k^ on this point, because aj^-^ becomes 
zero. 

Dispersion on a fc-line is proportional to a wave number for all direction in plane being perpendicular to the fc-line 
except for some singular points. From this results, fc-line is regarded as locus of a point where conduction and valence 
bands degenerate with fc-linear dispersion, i.e. the K (or K') point of the BZ in the graphene. This is a reason why 
we use a word " fc-line" . 



IV. DISCUSSION 

In this paper, we use a single-band tight binding model 
for explaining nature of networks which satisfy the con- 
ditions described in the section IIB. In the model, we 
assume that the hopping integral, is the same in any 
direction, in order to make a simple argument. But, this 
assumption is not necessarily required for networks hav- 
ing edge states. In a selected class of networks, we can 
prove the existence of edge states on them (See appendix 
A). _ 

Here, hypergraphite is defined to be a network which 
satisfies following two conditions. One is on the topology 
of a network given by hopping integrals and the other is 
on the electronic structure. 

The topology of network given by the hopping inte- 
grals is regarded as being constructed by the following 
method. 

1) Prepare an (TV — l)-dimensional AB bipartite net- 
work having a NBO whose amplitude is finite only 
on A-sites. The number of A-sites and that of B- 
sites in a unit cell are equal. 

2) Line up the copies of the network along z-direction. 
We define A-sites and B-sites on the i-th AB bipar- 
tite electronic network as the same way on the first 
network. 



3) Connect B-sites on the first network with A-sites 
on the second network by extra bonds. 

4) Connect B-sites on the i-th network with A-sites 
on the {i + l)-st network by the same way which is 
used in 3) (i = 1, 2, • • • , oo). 

Next, consider the s-TBM in the network, the elec- 
tronic structure constructed above fashion has to show 
following characteristics. 

1) Edge states appear in this network with a zigzag 
surface. 

2) The highest valence band and the lowest conduc- 
tion band degenerate at the Fermi energy and the 
Fermi surface becomes an {N — 2)-dimensional sur- 
face. 

3) The highest valence band and the lowest conduc- 
tion band possess fc-liner dispersion at points on 
the {N — 2)-dimensional Fermi surface. 

4) The bulk system is a zero-gap semiconductor. 

Here, note that the Fermi level is set on the center of 
the band structure, because we consider the case of a 
half-filled AB bipartite network. 

From the above definition on a network, the appear- 
ance of edge states on hypergraphite itself is natural. 
This is because a NBO which exists on a basic AB 
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bipartite network coincides with a completely localized 
edge states in a hypergraphite. But it is not clear 
whether hypergraphite become a zero-gap semiconduc- 
tor and whether it has bands with /c-linear dispersion 
at the Fermi energy. This important point is left to be 
solved in general. 

Here, we list possible scenarios how hypergraphite net- 
works and their characters are realized in real materials. 

We have shown that the three-dimensional three-fold 
coordinated network and the diamond structure are clas- 
sified as the hypergraphite. From the point of geometric 
topology, the cubic diamond and the cubic silicon possess 
the topology of the diamond structure. We also find that 
the 3-3 network is created by silicon in Q!-ThSi2. 

In these real materials, the electronic and geometric 
structures are determined by not only s-electrons but 
also p-electrons of carbon or silicon atoms. Actually, sp^ 
or sp^ hybridized orbitals form these interesting lattice 
structure and electronic properties. Hence, no electronic 
characteristic of hypergraphite is found around the Fermi 
energy. But, we notice the characteristics of the hyper- 
graphite in the lowest two bands. Examples are 1) gap 
closing at X-point in the lowest two sp^-hybridized bands 
of diamond or siliconEil and 2) gap closing at a point on 
the F-X line in the lowest two ^aap^-hybridized bands of 
CaSi2 in the a-ThSi2 structurecj. This is because the 
two bands are mainly composed of 2s or 3s electrons. 



V. SUMMARY 

We propose a class of A^-dimensional networks. These 
networks have the same characteristics of topology of 
hopping integrals and an electronic structure that the 
graphene has. A characteristic of topology of hopping in- 
tegrals is AB bipartite network. Particular properties of 
the electronic structure are 1) appearance of edge states 
under a zigzag surface and 2) existence of a fc-line with fc- 
lincar dispersion, and 3) being a zero-gap semiconductor. 
From the above reason, we regard the class of networks 
as an extended graphitic network and name them "hy- 
pergraphite" . 
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APPENDIX A: TRANSFER MATRIX METHOD 



In this appendix, we show that the edge states appear in a selected class of A^-dimensional networks which are 
composed of {N — l)-dimensional AB bipartite network. Naturally, we consider the TV-dimensional networks which 
satisfy the conditions of hypergraphite. 

We consider a semi-infinite slab model, on which the single-band tight- binding model is constructed. In this model, 
hopping integrals do not have to be unity. We call an i-th {N — l)-dimensional network A*, (i — 1,2, ■ ■ ■ ,n, ■ ■ ■ , oo) It 
is assumed that there exists a NBO on A^. We consider only AB bipartite networks, where 1) the number of A-sites 
and that of B-sites are the same and 2) the NBO has a finite amplitude on every A-site of A^ and all of B-sites are 
node for the NBO. Since all of A* has the same network as that of A^, all of A* has a NBO which is the copy of the 
first NBO. Then, we naturally define A-sites and B-sites on A' so that B-sites are always node for the NBO. A surface 
of the network is composed by only A-sites on A^. 

We also assume that the (n — l)-dimensional AB bipartite network is a lattice having a unit cell and is periodic 
in all (n — l)-directions. The numbers of A-sites and B-sites in the unit cell of A* are equal and represented by A^o- 
Then, we can obtain wave functions on A' in a Bloch form as, 

l^^(k)> -EE^^^-^'(k)exp(zk.rOc^_^.jO). (Al) 

3,1 r. 

Here k and represent an (N — l)-dimensional wave vector and a position vector for a unit cell on A^, respectively. 
Two labels, j and I, indicate sublattices {j = A or B) and a site of the j-sublattice in the cell {I = 1,- ■ ■ ,Nq), 
respectively. An operator, c]^, j i, creates an electron on a site indexed by j and / in a cell at on A*. 

We introduce an A^o-dimensional vector u^^ = (wij-_;(k)). Then we can rewrite the Schrodinger equation on the 
total A^-dimensional network in a form as, 

e u,.^ = T'"u,.B+T^-u,_i,B, (A2) 
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Here, matrices, T™ and T*^^, are functions of k. An (A*o x Nq) matrix represents the transfer from B-sites to 
A-sites in A*. T'^^ is another (A'o x Nq) matrix representing transfers from B-sites in A'~^ to A-sites in A*. T™^ is a 
matrix representing transfers from A-sites to B-sites in A* . 

In the first place, we discuss the simplest case, i.e. the unit cell in A* has only one A site and one B site. Then 
T™, T'^^ and Uij become a scalar, (i = 1, 2, • • • , n, • • • oo, j = A or B). We seek for a dumping wave with e = 0. 
Then we have the following equations with a transfer matrix, Tj, 

"^:r)^^(s)^(rr)(::;:)- 

Ui.A 



with another condition, ui^b = 0. Here, we assume that T'^^ and T™ are not equal to zero. This is natural because, 
a) if T*" is zero, there exists a completely localized edge state for a given k, and b) if T^^ = 0, there is no solution 
with ui^A 7^ and = 0. The problem is that on what conditions this set of equations has a dumping wave toward 
the z-direction. If we multiply two transfer matrices, we have, 



U'i+l,B \ _ f \ f Ui^B 



Ui+i.A J \ J V ^hA 



(A5) 



Then, we can easily see that a solution can be chosen to satisfy either Ui,A = or = 0. Because wi,_b = to 
satisfy the boundary condition of edge states, we have a solution, 

= ^ T'^u.^A. (A6) 

Thus, in region of reciprocal lattice space where k satisfies |7^«"/ye^| < have a dumping wave. The dumping 

factor D{k) is given by T'" /T" 



Next, we have to show that a wave vector satisfying |£'(k)| < 1 exists in a given network whose T*", T'^^ and 



'1,3 

are a scalar. But, we have assumed that a zero energy states exist on A^, which is given by ko satisfying T™(ko) = 0. 
Namely, this state is a localized eigen-state. Thus existence of a solution with D = is assumed. In usual networks 
with finite number of bonds for each site, T™ and T'^^ are analytic functions of fcj {j — 1, - ■ ■ , Nq). Hence, one can 
assume the continuity of T™ (and T'^^) as a function of kj. Then we have solutions with < \D\ < 1, which gives 
the degenerate edge states. 

The graphene and the 3-5 network shown in this paper correspond to this simplest case. In case of graphene, 
T'^^ = — 1 and = — 2cos-|. Hence, the dumping factor D becomes 2 cos I. While, in case of the 3-5 network, 

T"^ = -1 and T™ = -2(cos ^ + cos ^). Hence, the dumping factor D becomes D = 2(cos ^ + cos ^). 

So far, we have assumed that j^m scalars. In general, Ts and Uij could be a matrix and a 

vector, respectively. Then there is a case which T'^^ become a singular matrix. Hence, it is impossible to derive 



analogue of Eqs.(A3) and (A_4) from Eq.(|A2|). But, we may be able to obtain analogous determination equations for 
an A'^-dimensional network. 

In case when (yex-^-i g^ist, we see that an eigenvector of each A-sites on A^ is given by D^^^ui^a and Ui,B — 0. 
Thus, the eigenvector of each A-sites on A^ is represented by eigenvector of each A-sites on Ai_i and a dumping factor 
D as, 

Ui+i,A = -D Ui^A- (A7) 



Substituting Eq.([A7|) to Eq.(|A2|), we have, 

r"u,,A = -D T^^u,,^. (A8) 

The dumping factor, D, is determined by this equation. We can judge whether edge states emerge at an edge of the 
A^-dimensional network from it. When D — 0, the edge state is completely localized at the surface. When < \D\ < 1, 
the edge state becomes a dumping wave. When \D\ — 1, the state coincides with a bulk state. When 1 < \D\, the 
solution becomes unphysical. 

We apply this formulation to the 3-3 network. The equation becomes 

(l-t- exp(2ifca;) -exp(i(fca; + fcj^)) \ / ui^a \ ^_r,( ^ \ ( "1-4 
-1 -(1 + exp(2ifcy)) ) \ U2,A ) "^y-l Q )\ U2,A 

With this equation, we obtain a dumping factor for the 3-3 network. It is given hy D — Acoskx cos ky. 
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